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Abstract 

The Sasa-Satsuma equation is an integrable higher-order nonlinear Schrodinger equation. 
Higher-order and multicomponent generalisations of the nonlinear Schrodinger equation are 
important in various applications, e.g., in optics. One of these equations is the Sasa-Satsuma 
equation. We present the binary Darboux transformations for the Sasa-Satsuma equation and 
then construct its quasigrammians solutions by iterating its binary Darboux transformations. 
Periodic, one-soliton, two-solitons and breather solutions are given as explicit examples. 
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1 Introduction 

The celebrated nonlinear Schrodinger (NLS) equation 

iqt + ^qxx + \q\ 2 q = 0 (1.1) 

is considered to be one of the fundamental integrable equations admitting an n—soliton solution. It 
is proved integrable via the inverse scattering transform [28]. The NLS equation has applications in 
a wide variety of physical systems such as water waves [1,2,29], plasma physics [30] and nonlinear 
optics [11,12], This equation can be used to model short soliton pulses in optical fibres [16]. The 
basic phenomena are described by the nonlinear Schrodinger equation (1.1), but as the pulses get 
shorter various additional effects become important. In light of this fact, Kodama and Hasegawa 
[17, 18] developed a suitable higher-order NLS equation to take these additional effects. Their 
equation takes the form 

iqr + aiqxx + < 22 1 q\ 2 q + iS [fiiqxxx + h\q\ 2 qx + fizq (|<?| 2 ) x ] = 0 , ( 1 . 2 ) 

where the on, fii are real constants, <5 is a real spectral parameter and q is a complex valued function. 
The first three terms (Setting 5 = 0 and a 2 = 2oq = 1 ) form the standard nonlinear Schrodinger 
equation (1.1). In general, the Kodama-Hasegawa higher-order NLS equation (1.2) is not integrable 
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unless some restrictions are imposed on the real constants /% (i = 1, 2,3). With appropriate choices 
of these real constants, the inverse scattering transform can be applied to verify integrability of 
the resulting equation. It is known until now, along with the NLS equation (1.1) itself, there are 
four cases in which integrability can be proven via inverse scattering transform. These are the 
Chen-Lee-Liu [3] derivative NLS equation (f3\ : @2 : /?3 = 0 : 1 : 0), the Kaup-Newell [14] derivative 
NLS equation (/3i : fa : fa = 0 : 1 : 1), the Hirota [13] NLS equation (fa : fa : fa = 1 : 6 : 0) and 
the Sasa-Satsuma [24] NLS equation (fa : fa '■ fa = 1 : 6 : 3). 

Sasa and Satsuma [24] consider the case where aq = \ and a.^ = 1, that is 

iq.T + \<1XX + \q\ 2 q + iti [qxxx + 6\q\ 2 qx + 3 q (|g| 2 ) x ] = 0. (1.3) 

Sasa and Satsuma [24] introduce variable transformations 

»(*,*) = g(X,r) exp { =2 (x-T)j, t = T, x = X - T. (1.4) 

Then, setting 5 = 1, the equation (1.3) is reduced to a complex modified KdV-type equation 

u t + u xxx + 6] u\ 2 u x + 3 u (| u\ 2 ) x = 0 (1.5) 

which is an equivalent version of (1.3). The equation (1.5) is commonly known as the Sasa-Satsuma 
(SS) equation, and we will denote it as such from now on. The SS equation (1.5) is one of the 
integrable extensions of the NLS equation (1.1). The integrability of this equation has been widely 
studied with various methods such as the inverse scattering scheme [24], the Hirota’s bilinear 
approach [7,8] and the Darboux-like transformations [25]. 

In the present paper, we present a systematic approach to the construction of (1.5) by means of 
a standard binary Darboux transformation (BDT) and written in terms of quasideterminants [5,6]. 
Quasideterminants have various nice properties which play important roles in constructing exact 
solutions of integrable systems [9,10,19,22,27]. 

In this paper, we establish for the first time a standard BDT for the SS equation (1.5). Our 
solutions for the SS equation are written in terms of quasigrammians rather than determinants. It 
should be emphasised that these quasigrammian solutions arise naturally from the binary Darboux 
transformation we present here. Furthermore, we present periodic, one-soliton, two-solitons and 
breather solutions for the SS equation. 

This paper is organized as follows. In Section 1.1 below, we give a brief review on quasideter¬ 
minants. In Section 2, we construct a 3 x 2 eigenfunction and corresponding constant 2x2 square 
matrix for the eigenvalue problems of the SS equation (1.5) via a symmetry matrix. In Section 
3.2, we state a standard binary Darboux theorem for the Sasa-Satsuma system. In Section 3.3, 
we review the reduced binary Darboux transformations for the SS equation, which can be consid¬ 
ered as a dimensional reduction from (2 + 1) to (1 + 1) dimensions. In Sections 4, we present the 
quasigrammian solutions of the SS equation by using the binary Darboux transformation. Here, 
the quasigrammians are written in terms of solutions of linear eigenvalue problems. In Section 5, 
periodic, one-soliton, two-solitons and breather solutions of the SS equation are given for both zero 
and non-zero seed solutions. The conclusion is given in the final Section 6. 
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1.1 Quasideterminants 

In this short section we recall some of the key elementary properties of quasideterminants. The 
reader is referred to the original papers [5,6] for a more detailed and general treatment. 

Quasideterminants were introduced by Gelfand and Retakh in [5] as a natural generalisation of 
the determinant to matrices with noncommutative entries. Many equivalent definitions of quaside¬ 
terminants exist, one such being a recursive definition involving inverse minors. Let M = ( mij) be 
an n X n matrix with entries over an, in general non-commutative, ring 5^. Then the quasidetermi- 
nants of M for i. j = 1 ,,n are defined by 

\M\ij = mt, - r? (MV)' 1 ( 1 . 6 ) 

where rj is the row vector obtained from i th row of M with the j th element removed, c* is the column 
vector obtained from j th column of M with the i th element removed and MV is the (n — 1) x (n — 1) 
submatrix obtained by deleting the i th row and the j th column from M. Quasideterminants can 
also be denoted by boxing the entry about which the expansion is made 


I M\ij 





(1.7) 


If the entries in M happen to commute, then the quasideterminant \M\ij can be expressed as a 
ratio of determinants 


\M\ij 


(-1 ) i+j 


det M 
det Mv 


( 1 . 8 ) 


In this paper, we will consider only quasideterminants that are expanded about a term in the 

last column, most usually the last entry. For example considering a block matrix M = 

where A is an invertible (square) matrix over of arbitrary size and B, C are column and row 
vectors over 5^ of compatible lengths, respectively, and d E the quasideterminant of M is 
expanded about d is 


A B 
C d 


= d-CA~ l B. (1.9) 


Later we will use the following invariance of quasideterminants which follows immediately from 
their definition. Let a and /3 be invertible matrices of the same dimensions as A. Then 


A 

B 

C 

d 


aA/3 aB 


A 

B 

cp 

d 


C 

d 


( 1 . 10 ) 


2 The eigenvalue problems for the Sasa-Satsuma Equation 

The Lax pair [24] for the Sasa-Satsuma equation (1.5) is given by 

L = d x + JX + R (2.1) 

M = d t + 4JA 3 + 4RA 2 - 2QX + W, (2.2) 
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where J, R , Q and W are 3x3 matrices such that 


J = 


( 1 

0 

° \ 


0 

i 

° 

II 

cel 

V 0 

0 

-< / 




Q = 


\u\ 2 

u 2 

V'X \ 


u* 2 

\u\ 2 

i 

(2.3) 

K 

V*X 

-2 m 2 J 



and 


W = 


U U x — uu x 

0 

-4 u*\u\ 2 - u* 


uu x — u u x 
—4u\u\ 2 — u x 


ku\u\ 2 + U x , 
Au*\u\ 2 +u* x 

0 


(2.4) 


Here A is a spectral parameter and asterisk denotes complex conjugate. It can be seen that the 
potential matrix R in (2.3) has two symmetry properties [15,26]. One is that it is skew-Hermitian: 
R) = —R. The other one is that SRS -1 = R*. where 


S = S~ 1 = 



(2.5) 


Let (j) = {4>i,4 , 2,4 > 3) T be a vector eigenfunction for (2.1)-(2.2) for eigenvalue A. Using the second 
symmetry, it is easy to see that Scf>* = (</>?;, {/>*, (j)^) T is another eigenfunction, for eigenvalue —A*. 
Using these vector eigenfunctions we can define a 3 x 2 matrix eigenfunction 6 with 2x2 eigenvalue 
A 


satisfying 


<t> 1 4> 2 

fa <t>\ 

4>3 4>*3 


A = 


e x + JO A + R6 = o, 

n i A T/l A 3 i A 7“)/l A 2 


A 0 
0 -A* 


( 2 . 6 ) 


(2.7) 

( 2 . 8 ) 


3 Darboux transformations and Dimensional reductions 

3.1 Darboux transformation 

Let us consider the linear operators 


L = d x + ^2uid z y , M = d t + ^2vid l y , (3.1) 

Z =0 2—0 

where Ui,Vi € where 2^ is a ring, in general non-commutative. The standard approach to 
Darboux transformations [4,20,21] involves a gauge operator Gq = 9d y 9~ 1 , where 9 = 9(x,y,t ) is 
a solution to a linear system 


= M(<f>) = 0, 


(3.2) 
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where 4> is any eigenfunction of L and M, has the property of leaving the above linear problems 
invariant: 

m=m= o, 

where 0 = Gq(( j>) is an eigenfunction of L and M (3.1) and the linear operators L 
M = GqMGq 1 have the same forms as L and M : 

n n 

L = d x + ^2 Uid l y , M = d t + Vid l y . 

i=0 i —0 


3.2 Binary Darboux transformation 

Let Gq and G g be two standard Darboux transformations map two linear operators L and L onto 
a common linear operator L such that 


L = GqLGq 1 = GgLGj 1 . 

(3.5) 

Then one may define a binary Darboux transformation [21,23] B gg = G- 1 Gq such that L = 
B g qLBq 1 -. In order to define G g one needs an eigenfunction of L. This problem can be got around 

by using the formal adjoint operator constructed according to the rule (a9*)^ 
denotes the Hermitian conjugate of a. If L(r/>) = L(Q) = 0 and = L^(p) 

eigenfunction 6 as 

= {—dyf a\ where 
= 0, we derive the 

0 = 6>0(6», /0) _1 , 

(3.6) 

where the eigenfunction potential D is defined such that Q(0,p) y = p^9. We can 
binary Darboux transformation explicitly as 

then construct the 

B d:P = i-en(e, 

(3.7) 

with adjoint 


B~l = i-pn(9,p)- ta-V. 

(3.8) 

It should be pointed out that the above formulas obtained for Bg )P and B g the transformation 
makes sense for any m x n matrices 6 and p such that L{6 ) = L^(p) = 0, not just square m x m 
matrices. Here, only the eigenfunction potential Q(9,p) must be an invertible square matrix. 

Suppose cip] = <f> is a general eigenfunction of the operators Lp] = L, Mpj = M and ^[i] = V 7 a 
general eigenfunction of the adjoint Lax operators = L\ = M\ where L and M are given 

by (3.1). We then define the binary Darboux transformations of the eigenfunctions 4> and V 7 as 

0[2] = B e, P (^[i]) = <%] - 0[i]^ (%],P[i]) 1 ^ (<£[ 1 ],P[ 1 ]) ; 
V’p] = B e,l ~ p[i] n i 9 [i]^[i])~^ n ^ 

(3.9) 

(3.10) 

with 


0[2] = 0[2]|<£-+0 2 i P[2] = "0[2] ly> — >-p2 ' 

(3.11) 
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(3.3) 

= GqLGq 1 and 

(3.4) 


After n ^ 1 iterations, the nth binary Darboux transformation is given by 

4 > [n+ 1] (^[n]) 4*[.n] ^[n]^ 1 (^[n]?P[n]) ^ i^[n\'> P[n\) 5 (3.12) 

0[n+1] = (0[n]) = 0[n] - P[ra]^ (#[n], P[n]) f ^ (#[n] , 0[n]) f > (3.13) 

with 

®[n\ — 0[n] !(/)—>■ 0 n > P[n] V^n] If/’^-Pn > (3.14) 

where L(9i ) = M(9i) = 0 and 10 (p*) = M^(pi) = 0. Using the notation 0 = ( 60 ,..., 9 n ) and 
P = (pi, ■ ■ ■, Pn) i we have 


and 


>[n+l] 


51(0, P) 
0 


51(0, P) 


V’[n+1] — 


51(0, P)t 51(0,0) f 

p ny 


^(0[n+l], 0[n+l]) 


51(0, P) 

51(0, P) 

51(0,0) 

51(0,0) 


(3.15) 


(3.16) 


3.3 Dimensional reductions of the binary Darboux transformation 

Here, we describe a reduction of the binary Darboux transformation from (2 + 1) to (1 + 1) dimen¬ 
sions. We choose to eliminate the y-dependence by employing a ‘separation of variables’ technique. 
The reader is referred to the paper [23] for a more detailed treatment. We make the ansatze 

0 = 4> r (x,t)e Xy , 9 = 9 r (x, t)e Ay , (3-17) 

0 = 0 r (x, t)e ra , p = p r (x,t)e Uy , (3.18) 

where A, p are constant scalars and A, n are N x N constant matrices and the superscript r denotes 
reduced functions, independent of y. Hence in the dimensional reduction we obtain <9* (0) = A*0 
and dy (9) = 9A 1 and so the operator L in (3.1) becomes 

n 

V = Qr + X^A 4 , (3.19) 

i =0 

where 9 r is a matrix eigenfunction of L r such that L r ( 9 r ) = 0, with A replaced by the matrix A, 
that is, 

n 

«£ + X u i 9 r A i = 0. (3.20) 

i=0 

It follows that the y— dependence of the potential 51 can also be made explicit by setting 51 (9,p) = 
e nty 51 r (9 r , p r ) e Ay and 51(0, p) = e( nt+A/ )^51 r (0 r ,/3 r ). Then the dimensionally reduced binary 
Darboux transformations are written as 

Bgr tP r = i-e r n r (9 r ,p r )- 1 (^ + ^i)~ 1 p r \ (3.21) 
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B ~t — I — 

1J e r ,p r ~ 1 


P r n r (d r ,p r )~^ (A t + / i/) e r \ 

where fl r is an algebraic potential satisfying the following conditions 

u^a r (0 r ,p r ) + n r (0 r ,p r )A = p r ^e r , 

(nt + a/) Q r (pf, f f) = P r V, 

n r (9 r ,ip r ) (a + /At) = i’ r] e r . 

The transformed operators 

rr _ jD rr o—l 

lj [n+ 1] — ^O r ,P r - Ll [n}^e r ,pri 

M [n+l] = 

have generic eigenfunctions and adjoint eigenfunctions 


£ 


n+l] 


-t 


>P[: 


n > U n / > 


V’E.+l] = d»l - 4.1^ («w . P\n l) _ «’ («W' d»l) 


with 


^[n] - 0[n] 1 ) P[n] “ Vf n ] I V^PS.' 

From now on, for notational simplicity, we omit the superscript r. 


(3.22) 

(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 


4 Quasi-Grammian solutions of the Sasa-Satsuma equation 

In this section we determine the effect of the BDT Bg )P = I — 00(0, p) -1 (III + A/) 1 pi on the 
operator L = d x + J A + R given by (2.1) with 9 an eigenfunction of L and p an eigenfunction of 
Ll. Corresponding results hold for the operator M given by (2.2) and its corresponding adjoint 
M 1. In the operator L, the matrix coefficients J and R are both skew-Hermitian. This leads us to 
observe the relation L + = 0 in which we can let p = 9 and ip = <p. 

The operator L is transformed to a new operator L such that 

L = BgLB~\ (4.1) 

where Bg = I - 00(0, 0)" 1 (ill + a/) - 1 0t. 

We find that 


R = R + 


.7,0fi(0,0) _1 0l 


(4.2) 


For notational convenience, we introduce a 3 x 3 matrix P = (pij) such that R = [J, P], and hence 



p 11 

P 12 

U 

P 21 

P22 

u* 

u* 

U 

P33 


(4.3) 
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From (4.2), since R = [J,P], it follows that 


P = p-en(e,e)- 1 e\ ( 4 . 4 ) 

in which the relation (3.23) is now 

0(0, 6)A - A t f7(0, 6) = 9^9 (4.5) 


with II = —A, where the eigenfunction 9 and the diagonal constant matrix A are given by (2.6). 

Let = P, P[ 2 ] = P, 9[ i] = 9i = 9 and Ai = A so that Ai = diag (Ai, — A^). Then, the 
solution (4.4) can be rewritten as 

P[2\ = P[i] ~ <9[i]^(%],<9[i]) _1 ^[i]- (4.6) 

After n repeated applications of the reduced binary Darboux transformation Bg, we have 

P[n+1] = P[n] ~ 0 N^(^[n], ^[n]) _1 ^[ n ] (4-7) 

with = 0[ n ] | ( j,^e n , where the general eigenfunction 0[„] is given by (3.28) as 

4 > [n+ 1] ^[n] ^[n]^(^[ra] > ^[n]) ^(^[n] > ^[ra]) • (4.8) 

Let 0i,..., 4> n be a particular set of eigenfunctions of the linear operators L, M given by (2.1) — (2.2), 
and define 0 = (9\, ..., 9 n ) for the 3x2 matrices 9i (i = 1 ,..., n) such that 


( foi-2 03.-1 \ 

= j 03i— 1 03i_2 I • (4-9) 

V ^3* 03 i / 


We express P[„+i] and 0[ ri+1 ] in quasi-Grammian forms as 


P[n+ 1] 



0(0,0) 


0t 


p + 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 


0[n+1] 


0 ( 0 , 0 ) 

0 


0 ( 0 , 0 ) 

0 


(4.10) 


In order to express the quasi-Granunian solution Pr n+1 i explicitly in terms of the variable u, we 
define the 3 x 2 n matrix 0 as 


0 = 



(4.11) 


where dq, 4> 2 and 4> 3 denote the row vectors ( 0 i, 02,, 03n-2, 0 L-i)> ( 02 , 0 *, ■ • •, 03n-i, 03 n _ 2 ) 
and ( 03 , 03,..., 03 n , 03 n ) respectively. Thus, we obtain 









p 


71+1] 


= P + 


12(0,0) 

$1 fo 


n(0,0) 




«(©,©) <*4 


$i 


o 


0 ( 0 , 0 ) 

$2 [0 

0(0,0) 4*{ 

cj) 3 


0 ( 0 , 0 ) 

$2 [F 

0 ( 0 , 0 ) 

IF 


0(0,0) $1 
4>2 |F 

0 ( 0 , 0 ) 

$3 [F 


(4.12) 


By substituting (4.3) into (4.12), we have quasi-Grammian expressions for u and u*, namely 


U[n+ 1] = U + 2i 

= u + 2i 


0 ( 0 , 0 ) ^ 


$i 


0 


0 ( 0 , 0 ) 

$3 IF 


U I 


71+1] 


= w* + 2i 


= ii* + 2i 


0 ( 0 , 0 ) 

4> 3 IF 


0 ( 0 , 0 ) $1 

$2 |0 


(4.13) 

(4.14) 

(4.15) 

(4.16) 


These form a quasi-Grammian solution u of the Sasa-Satsuma equation (1.5) and its complex 
conjugate, However, it is necessary to show that these four expressions are consistent. That is, that 
the expressions on the right hand sides of (4.13)—(4.14) and (4.15)-(4.16) are equal and that the 
pairs are indeed complex conjugates. The proof of this is given in Section 4.1. 


4.1 Proof of consistency 

In the expressions (4.13)—(4.16), the potential 12(0, 0) is a 2n x 2 n matrix satisfying the relation 


12(0, 0)A — A t P(0, 0) = 0 t 0, (4.17) 

where A is 2n x 2 n constant matrix such that A = diag (Ai,..., A n ). Solving this relation for 
12(0,0), we obtain the explicit expression 


/fi(0!,0l) ^(02,01) 

n(d 1 ,e 2 ) n(e 2 ,e 2 ) 


12 (#n, #l)\ 

^(9 n ,9 2 ) 


yi(0 i,Q n ) Q(9 2 ,9 n ) ... Q(9 n ,0 n )J 


(4.18) 


where 12 (#j, Oj) is 2 x 2 potential satisfying the relation 

n(e u d j )A i -A]n(e il 9 J ) = 0]e i , 


(4.19) 
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where = diag (A*,, —Ajjl) and i. j. k e {1,2,..., n}. It follows from this relation that the potential 
O can be written explicitly as 


where Fij = Fij(x, t, A*, Aj) and G^ = Gij(x, t, A ? ;, Aj) are the scalar functions 

1 


(03i-2<^3i-2 + foi-l&tj-l + <t>Zi4>lj) , 


- a,- y 

Gjj = —— (03j_203i-l + 03i—103jr —2 + 4>3i4>3j) ■ 

Aj ~r Aj 


(4.20) 

(4.21) 

(4.22) 


Here we observe that F^ and Gjj are such that F*- = — Fji and Gij = Gjj, for i, j = 1,..., n. Then 
the 2x2 potentials 0(6{, 6j) satisfy the symmetry condition 


£l(9i,0j) + £1(0j, 0^ = 0, 

and the 2 n x 2 n matrix potential 0(0, 0), as given by (4.18), is skew-adjoint, 

0 ( 0 , 0 ) + 0 ( 0 , 0 )* = 0 . 


Now it is readily seen that 


0(0,0) d>3 




0(0,0)t 

4-3 [0 


0 ( 0 , 0 ) 

*1 

$ 3 

0 


(4.23) 

(4.24) 

(4.25) 

(4.26) 


using (4.24), and so (4.13), (4.15), and similarly (4.14), (4.16), are indeed complex conjugate. 
It remains to prove that (4.13) and (4.14) are consistent, i.e. that 


0(0,0) 

4 


0(0,0) 

<4 


0 


cj>3 

0 


(4.27) 


Note first that this each side in this equation represents a scalar and so the right hand side can 
also be written as 


0(0,0) 

*5 

(J> 3 

0 


0 ( 0 , 0 ) T 

€ 

^*2 

0 


(4.28) 


where T denotes the matrix transpose. 

Now let a be the 2 n X 2 n permutation matrix 


/0 1 ••• 0 0\ 

1 0 ••• 0 0 


0 0 ••• 0 1 

\0 0 10/ 


(4.29) 
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Pre(post)-multiplying any matrix with 2 n (rows) columns by a has the effect of interchanging its 
2ith and (2 i + l)th (rows) columns for i = 1,..., n. Hence 


4qa = 4>£, 4* 2 a = <h 3 a = $3 

1 

aO(0,0)a = 0(0, 0) T . 

Now using the quasideterminant invariance (1.10), we get 


0 ( 0 , 0 ) 



aO(0,0)a a4>3 

$1 

0 


4>i a 

0 




H( 0 , 0 ) T <I >3 





4)2 0 



This completes the proof. 


(4.30) 

(4.31) 


5 Particular solutions 


In order to construct particular solutions for the Sasa-Satsuma equation (1.5), we consider the 
quasi-Grammian solution given by (4.13) 


u [n+ 1] 


u + 2i 


0(0,0) 


4>i 

0 


(5.1) 


where 4>i and $3 denote the first and third rows repectively of a 3 x 2 n matrix eigenfunction 0, 
and the potential 0(0,0) is the 2 n x 2 n matrix with entries defined by (4.18) and (4.20)-(4.22). 

Let us consider the spectral problem L(0) = M(^) = 0 with eigenvalue A j (j = 1,..., n ), where 
4> = <t> 3 j-i, <t> 3 j) and L and M are given by (2.1)-(2.2) so that 


<t>x + J <^Aj + R<p — 0, 
(p t + 4J(pX!j + 4B4X 2 3 - 2Q(pX :j + W(t> = 0, 

where J, R, Q and W are given by (2.3)-(2.4). 


(5.2) 

(5.3) 


Case (n=l) 

In this case 4> = (<^> 1 , cf) 2 , (t>z) T is a solution of the spectral problem L(<)>) = M(cj)) = 0 with eigenvalue 
Ai. Thus, from (5.1), we derive the following explicit solution 


u [2] ~- 

where 

Tn = 

G n = 



Tn 

-c\ 1 

</>3 

u + 2i 

Gn 

Tn 

<t> 3 


<t>i 

to* 

0 


Ai ]_ (l^il 2 + \H 2 + Ifol 2 ) , 

( 201^2 + 03 ) • 


(5.4) 


(5.5) 

(5.6) 
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5.1 Solutions for the vacuum 

For u = 0, the above eigenvalue problems (5.2)-(5.3) transform into the first-order linear system 


4*x T J 0Aj — 0, 

<t>t + 4 J (j)Xj = 0, 

which has solution (j) = (</>i, </>i, ipi) T , where 

</,j = e ~ iX ^ x+4X ^\ ipj = ( iX A'-' AX V)- 


(5.7) 

(5.8) 


(5.9) 


Case (n=l) Then, the solution (5.4) becomes 


up] = 2 i 


f -9* i’l 

g f i>i 


4>i <j>\ [o 


which can be written as 


where 


U[2] = — 2 i 


f = 


. / Ai + \ g + g 


Ai-A \) p + \g | 2 ’ 

(2|0 1 | 2 + |^i| 2 ) 


Ai-AJ 

1 


9 — (^l + ^i) ' 


Here the eigenfunction cj) = in which 

</>! = e -*Ai(x+4A \t), ^ = e iAi(x+4A ft). 

By setting Ai = £ + irj (£, rj £ M and £, g / 0), we have 

£(<7 + 9* 


U[2] 2 >itr 2 + i/y i 2 ) 


2V 


in which 


/ - + 


5 = 


_1_ ( 2 e 2 ( Q -*^) + g -2(a-i/3)\ 

+ in) V J 


(5.10) 


(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 

2K + ! „)W ■- r (5 ' 17) 

where a = g [x + 4[3£ 2 — r/ 2 ]t) and /? = £ (x + 4[£ 2 — 3r/ 2 ]f). Then we obtain a breather solution 
£ [3cosh(2a) + sinh(2a)] cos(2/3) — r] [cosh(2a) + 3sinh(2a)] sin(2/3) 


^ ^ £ 2 [3 cosh(2a) + sinh(2a)] 2 + 8r/ 2 sin 2 (2/3) 

The solution (5.18) is plotted in the Figure 1. 


(5.18) 
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Figure 1: Breather solution up] of the SS equation (1.5) with the choice of parameters £ = 0.4, rj = 
0.5. 


5.2 Solutions for non-zero seeds 

For u / 0, u = k is a solution of the Sasa-Satsuma equation (1.5), where we let Ho be a 
real constant. We use this as a seed solution for application of binary Darboux transformations. 
Substituting u = k into the linear system (5.2)-(5.3) and then solving for the eigenfunction 0 = 
(</) 3j _2,03j~i,03i) T , we obtain 


03.7-2 = -X 


037-1 - 


-cie _iAj '( x+4A ? t ) + C2e V A i +2fc2 [ a:+4 ( A ? _fc2 ) t ] + c 3 e“V^ +^[ x + 4 { x h k2 ) t \^ (5 19 ) 
1 c ie - iX ^ x+4X i t ) + c 2 eV A ^ +2fc2 [ a: + 4 ( A ? _fc2 ) t ] + ^ ( 5 20 ) 


031 = A0e 


i\/ A J+ 2 fc 2 [x+ 4 ( A ?-fc 2 )i] + A ^ e -V A2+2A;2 [ a;+4 ( A ?" fc2 ) t ]. 


(5.21) 


where A] = ^Aj + y^A 2 + /c 2 ^, A 2 = ^Aj — y^A 2 + fc 2 ^ and ci, c 2 , c 3 are arbitrary constants. 


Case (n=l) The solution (5.4) can be written as 


u \2] = k — 2i 


. (0105 + 0503) An + 0103^! - 0505Gn 


^i 2 i + |Gii| 2 

where An and Gn are given by (5.5)-(5.6). Substituting 0i,0 2 and 0 3 into (5.22), we obtain 


(5.22) 


where 


.. , ,,, (I0i| 2 - I0a| 2 ) (0105 - 0 2 0s) 

Tirol — rt t 4U-ft, 

(1011 2 - 102 | 2 ) 2 + 2 10105 - 05031 


01 = — ^ Cl X*- 4/A ) + C2 e iD l x - 4 ( k2+ ^) t \ + C3 e- iD M( fe2 +M 2 )*], (5.24) 

0 2 = 0i + ciH( a; - 4 0 !f ), (5.25) 

03 = -^-iD)^- 4 ^^ 2 )*] _ ^( /x + 5 j jD ) e -^M(fe 2 H 2 )t] ) (5.26) 

ft ft 

in which A = y / 2/c 2 — /z 2 . Here, for simplicity, we have chosen Ai E zM such that Ai = i/i. 


(5.23) 
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Figure 2 : Periodic solution uy 2 \ of the SS equation (1.5) with the choice of parameters k = 0.5, fi = 

0 . 1 . 


Case D 2 (n) = 2 k 2 — /i 2 > 0 

For the case D 2 = 2 k 2 — /i 2 > 0, (5.23) can be written as 

c 2 Ke 2ia + c* 2 K*e~ 2ia + 2/r|c| 2 - 2kS{ce ia + c*e~ ia )e~P 

U [2] -k- 2 kfi ^ c 2 Ke 2ia + c *2 K * e ~2 ia) + 4fc 2| c |2 + 2 fc 2 1S'| 2 6“ 2 /5 ’ 

where a = H[x — 4(/c 2 + /r 2 )t], /3 = n(x — 4/i 2 t), K = /i — iD, S = ^(|c 2 | 2 — |C 3 1 2 )and c = 
c\c 2 + C 1 C 3 / 0 . 

By choosing c 2 = C 3 so that 5 = 0, we obtain a periodic solution 


(5.27) 


ti[ 2 i — —k T 2- 


fc(2fc 2 — /r 2 ) 




^cos( 2 a) + \j2k 2 — /r 2 sin( 2 a) 


+ 2fc 2 


(5.28) 


This solution is plotted in the Figure 2. If we choose c\ = c 2 — C 3 then c = |c 2 | 2 — | C 3 1 2 and 
S = iicD. The solution (5.27) can be written in the following form 


u [2] = k — 2 k/jL- 


/rcos(2a) + D sin(2a) + n — 4iH[cosh/3 — sinh/3] cos a 


/u[/ucos(2a) + Z4sin(2a)] + 2 k 2 + 4H 2 [cosh(2/3) — sinh(2/3)] 
This solution is plotted in the Figure 3. 


(5.29) 


Case D 2 (n) = 2k 2 — n 2 < 0 

For the case D 2 = 2k 2 — /i 2 < 0, (5.23) can be written as 

_ 1 ji{e 2 ^+ e~ 2 ^)-E{e 2 ^-e~ 2 ^) + 2^ + 4:Eh^ 1 h 2 (e^+ e~^)e~^ 

U[2] M ^ 2 (e 2 ^ + e- 2 ^) - ^£(e 2 ^ - e” 2 ^) + 4A ; 2 - 8 E 2 h^ 2 h 2 e~ 2 P 


(5.30) 


where E = \J[i 2 — 2k 2 , 7 = E[x — 4 (k 2 + /x 2 )t], /3 = — 4^z 2 t) and /ii = cic 2 + c*C 2 7 ^ 0, 

h 2 = c 2 C 3 — c 2 C 3 such that h\ E M and h 2 E iM. 
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Figure 3: Solution |u [ 2 ]| 2 of the SS equation (1.5) with the choice of parameters k = 0.7,/r = 0.3. 



Figure 4: One-soliton solution |tt [ 2 ]| 2 of the SS equation (1.5) with the choice of parameters 
k = 0.1,/i = 0.5. 


By choosing C 2 = C 3 so that /12 = 0, we obtain 


u [2\ 


-k- 2 


k([i 2 — 2k 2 ) 

fj, 2 cosh(27) — f.fj, 2 — 2 k 2 sinh(27) + 2 k 2 


(5.31) 


This is one-soliton solution. This solution is plotted in the Figure 4. If we choose C 3 
h'2 = ih\. Then the solution (5.30) can be written as 

E(l — 4e -2/3 ) + 4i/rcosli7e _/3 

^ ' ' fj , 2 cosh(27) — [jlE sinh(27) + 2 k 2 + 4 E 2 e~ 2 P 

This is two-soliton solution. This solution is plotted in the Figure 5. 


—ic\ so that 


(5.32) 
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Figure 5: Two-soliton solution |u[ 2 ]| 2 of the SS equation (1.5) with the choice of parameters 
k = 1, fi = 1.5. 

6 Conclusion 

In this paper, we have presented a standard binary Darboux transformation for the SS equation 
(1.5) and using this we have constructed a wide family of solutions in quasigrammian form. These 
quasigrammians are expressed in terms of solutions of the linear partial differential equations given 
by (2.7)-(2.8). Moreover, periodic, one-soliton, two-solitons and breather solutions for zero and 
non-zero seeds have been given as particular examples for the SS equation. Examples of these 
solutions are plotted in the Figures 1-5 for particular choices of parameters. 

One should notice that we have chosen u = k, where k is a real constant, as a seed solution 
of the SS equation. This is the simplest non-zero seed. However, one might also choose the seed 
u = fe e - ia ( 3;+ [ a2 - 6fc2 ] t ), where a,k G R, to construct various rich explicit solutions including those 
we present here. Furthermore, in the present paper, we only consider the case n = 1 for constructing 
explicit solutions. It can be obtained more explicit solutions by considering other cases such as 
n = 2. Finally, it should be pointed out that the binary Darboux transformation technique is a 
universal instrument that allows us to construct exact solutions for other integrable systems. 
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